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Abstract 

A theory for constructing integrable couplings of soliton equations is developed 
by using various perturbations around solutions of perturbed soliton equations being 
analytic with respect to a small perturbation parameter. Multi-scale perturbations 
can be taken and thus higher dimensional integrable couplings can be presented. The 
theory is applied to the KdV soliton hierarchy. Infinitely many integrable couplings are 
constructed for each soliton equation in the KdV hierarchy, which contain integrable 
couplings possessing quadruple Hamiltonian formulations and two classes of heredi- 
tary recursion operators, and integrable couplings possessing local 2 + 1 dimensional 
bi-Hamiltonian formulations and consequent 2 + 1 dimensional hereditary recursion 
operators. 



Running title: Integrable couplings by perturbations I. 



1991 Mathematics Subject Classification: 58F07, 35Q53. 

Key words and phrases: Integrable coupling, Multi-scale perturbation, Hereditary recursion 
operator, Lax representation, Bi-Hamiltonian formulation, KdV hierarchy. 



1 



1 Introduction 



Integrable couplings are a quite new interesting aspect in the field of soliton theory Q. 
It originates from an investigation on centerless Virasoro symmetry algebras of integrable 
systems or soliton equations. The Abelian parts of those Virasoro symmetry algebras 
correspond to isospectral flows from higher order Lax pairs and the non-Abelian parts, to 
non-isospectral flows from non-isospectral Lax pairs g, ||. If we make a given system of 
soliton equations and each time part of Lax pairs of its hierarchy to be the first component 
and the second component of a new system respectively, then such a new system will keep 
the same structure of Virasoro symmetry algebras as the old one. Therefore this can lead 
to a hierarchy of integrable couplings for the original system. 

Mathematically, the problem of integrable couplings may be expressed as: for a given 
integrable system of evolution equations ut = K(u), how can we construct a non-trivial 
system of evolution equations which is still integrable and includes ut = K(u) as a sub- 
system? 

Therefore, up to a permutation (note that we can put some components of ut = K{u) 
seperately), we actually want to construct a new bigger integrable system as follows 

Ut = ^ U \ (1.1) 
v t = S(u,v), 

which should satisfy the non-triviality condition dS/d[u] ^ 0. Here [u] denotes a vector 
consisting of all derivatives of u with respect to a space variable. For example, we have [u] = 
(u,u x ,u xx , ■ • • ) in the case of 1 + 1 dimensions. The non-triviality condition guarantees 
that trivial diagonal systems with S(u, v) = cK(v) are excluded, where c is an arbitrary 
constant. 

There are two facts which have a direct relation to the study of integrable couplings. 
First, all possible methods for constructing integrable couplings will tell us how to extend 
integrable systems, from small to large and from simple to complicated, and/or how to 
hunt for new integrable systems, which are probably difficult to find in other ways. The 
corresponding theories may also provide useful information for completely classifying in- 
tegrable systems in whatever dimensions. Secondly, the symmetry problem of integrable 
systems can be viewed as a special case of integrable couplings. Strictly speaking, if a 
system of evolution equations ut = K(u) is integrable, then a new system (called a per- 
turbation system) consisting of the original system and its linearized system 

u t = K(u), 
v t = K'(u)[v}, 
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must be still integrable Q. The second part of the above new system is exactly the system 
that all symmetries need to satisfy, but new system itself is a special integrable coupling 
of the original system u% = K(u). Generally, the search for the approximate solutions 
un = ^2iLo £%r \i^ N > 1, oi physical interest to a given system u t = K(u) can be cast into 
a study of the general standard perturbation systems r\i t = \ 9 K q^ n ^ > < i < N . These 
perturbation systems were proved to form integrable couplings of the original system 
Ut = K(u) [f|, H|, the simplest case of which is the above system associated with the 
symmetry problem. This fact is also a main motivation for us to consider the problem of 
integrable couplings. 

However the standard perturbation systems above are just special examples of inte- 
grable couplings. They keep the spatial dimensions of given integrable systems invariant 
and only the perturbations around solutions of unperturbed integrable systems have been 
considered. We already know || that it is possible to extend the standard perturbation sys- 
tems and to change the spatial dimensions, in order to make more examples of integrable 
couplings. The question we want to ask here is how to do generally, or what related theory 
we can develop. In this paper, we would like to provide our partial answer to this extensive 
question, by establishing a theory on the multi-scale perturbation systems of perturbed 
integrable systems. An approach for extending the standard perturbation systems and for 
enlarging the spatial dimensions by perturbations will be proposed. 

Let us now introduce our basic notation and conception, some notation of which comes 
from Refs. |, g |, |]. Let M = M{u) be a suitable manifold possessing a manifold variable 
u, which is assumed to be a column vector of q functions of t G II and x G R p with t playing 
the role of time and x representing position in space. We are concerned with coupling 
systems by perturbations and thus need to introduce another bigger suitable manifold 
Mjy = ^n(vn) possessing a manifold variable r)N = (rjQ,r]J, ■ ■ ■ ,rj N ) T , N > 0, where r/j, 
< i < N, are also assumed to be column function vectors of the same dimension as u and 
T means the transpose of matrices. Assume that T(M),T(Mn) denote the tangent bundles 
on M and M N , T*(M),T*(M N ) denote the cotangent bundles on M and M N , and C°°(M), 
C°°(Miy) denote the spaces of smooth functionals on M and Mn, respectively. Moreover 
let TJ(M) be the s-times co- and r-times contravariant tensor bundle and TJ| U (M), the 
space of s-times co- and r-times contravariant tensors at u € M. We use X{u) (not X\ u ) 
to denote a tensor of I £ TJ(M) at u € M but sometimes we omit the point u for 
brevity if there is no confusion. Note that four linear operators <3? : T{M) — > T(M), 
* : T*(M) -> T*(M), J : T*(M) -> T(M), 9 : T(M) -» T*(M) may be identified with 
the second-degree tensor fields T$ G T^(M), G T^M), Tj G T 2 (M), T e G T 2 °(M) by 
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the following relations 

T^(u)(a(u),K(u)) =< a(u),$(u)K(u) >, a G T*(M), K G T(M), 
T^(u){a{u),K{u)) =< V(u)a(u),K(u) >, a G T*(M), K G T(M), 
rj(«)(a(«),/3(«)) =< a{u), J(u)P(u) >, a,/3 G T*(M), 
Tq(u)(K(m),S(u)) =< 6(m)K(m),5(u) >, K,S G T(M), 

where < ■, • > denotes the duality between cotangent vectors and tangent vectors. 

Of fundamental importance is the conception of the Gateaux derivative, which pro- 
vides a tool to handle various tensor fields. For a tensor field X G TJ(M), its Gateaux 
derivative at a direction Y G T(M) is defined by 

X> ( u)[Y(u)\ = aY( " + er( " >) . (1.2) 

0£ £ = 

For those operators between the tangent bundle and the cotangent bundle, their Gateaux 
derivatives may be given similarly or by means of their tensor fields. The commutator of 
two vector fields K,S G T(M) and the adjoint map adx ■ T{M) — ► T(M) are commonly 
defined by 

[K,S]{u) = K'(u)[S(u)} - S'(u)[K(u)}, &d K S=[K,S\. (1.3) 

Note that there are some authors who use the other commutator 

[K,S}(u) = S'(u)[K( U )}-K'(u)[S(u)}. 

It doesn't matter of course but each type has many proponents and hence one must be 
careful of plus and minus signs in reading various sources. 

The conjugate operators of operators between the tangent bundle and the cotangent 
bundle are determined in terms of the duality between cotangent vectors and tangent 
vectors. For instance, the conjugate operator & : T*(M) — > T*(M) of an operator $ : 
T(M) -» T(M) is established by 

< &{u)a{u),K(u) >=< a{u),<f>{u)K(u) >, a G T*(M), K G T{M). 

If an operator J : T*(M) -> T(M) (or 6 : T(M) -> T*(M)) plus its conjugate operator is 
equal to zero, then it is called to be skew-symmetric. 

Definition 1.1 For a functional H G C°°(M), ite variational derivative ^ G T*(M) is 
defined by 

fiff 

< —(u),K(u) >= H'(u)[K(u)}, K G T{M). 
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If for 7 G T*(M) there exists a functional H G C°°(M) so that 
fiff 

— = 7, i.e., tf'(u) =< 7 (u),if(u) >, if G T(M), 

£/ien 7 € T*(M) is called a gradient field with a potential H. 

A cotangent vector field 7 G T*(M) is a gradient field if and only if 
(d 7 )(u)(K(u),S(u)) 

:= < 7'(«)[if («)], 5(u) > - < 7'(n)[,S(n)], if(u) >= 0, K,S £ T(M). (1.4) 
If 7 G T*(M) is gradient, then its potential H is given by 

H(u) = / < j(Xu),u > dX. 
Jo 

Definition 1.2 For a linear operator <3? : T(M) — > T(M) and a vector field K G T(M), 
the Lie derivative Lk& '■ T(M) — > T(M) of <I> with respect to K is defined by 

{L K <&){y)S{u) = $'(u)[K(u)]S(u) - K'(u)[<S>(u)S(u)} + $(u)K'(u)[S(u)}, S G T(M). 

(1.5) 

An equivalent form of the Lie derivative is 

(L K <Z>)(u)S(u) = *(u)[K{u), S(u)} - [K(u),<Z>(u)S(u)}, (1.6) 
where $ : T(M) T(M), K,S G T(M), and the commutator [•, •] is defined by (O). 



Definition 1.3 A linear operator $ : T(M) — > T(M) is called a recursion operator of 
Ut = K(u), K G T(M), if for all S G T(M) and u e M, we have 

^^S(u) + &{u)[K(u)}S(u) - K'(u)[^(u)S(u)} + $(u)K'(u)[S(u)} = 0. (1.7) 



Obviously a recursion operator $ : T(M) — > T(M) of a system -uj = K(u), K G T(M), 
transforms a symmetry into another symmetry of the same system ut = K(u). Therefore 
it is very useful in constructing the corresponding symmetry algebra of a given system 
and its existence is regarded as an important characterizing property for integrability of 
the system under study. 
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Definition 1.4 A linear operator $ : T(M) — > T(M) is called a hereditary operator or to 
be hereditary if the following equality holds 

&(u)Mu)K(u)]S(u) - <f>(u)$'( U )[K(u)]S(u) 

-&(u)mu)S(u)]K(u) + ${u)&{u)[S{u)}K(u) =0 1 ' j 

for all vector fields K,S £ T(M). 



For a linear operator <E> : T(M) — ► T(M), The above equality ( |1.8| ) can be replaced 
with either of the following equalities: 

= $(«)(L K $)(u), K G T(M), 
$ 2 (n)[i^(u),5(n)] + [$(u)K(u),$(u)S(u)] 

-<f>{u){[K(u),<S>{u)S(u)} + [$(u)iT(u),S(u)]} = 0, K,S G T(M). 



It follows directly from ( |1.6| ) that these two equalities are equivalent to each other. We point 
out that hereditary operators have two remarkable properties. First, if $ : T{M) — » T{M) 
is hereditary and = 0, K G T(M), then we have [$ m if, $ n iT] =0, m, n > (see, 

for example, || 10, |TT|] ) . Therefore, when a system u t = K(u), K G T(M), possesses a 
time-independent hereditary recursion operator $ : T(M) — » T(M), a hierarchy of vector 
fields <& n i^, n > 0, are all symmetries and commute with each other. Secondly, if the 
conjugate operator ^ = & of a hereditary operator $ : T(M) — > T(M) maps a gradient 
field 7 G T*(M) into another gradient field, then ^"7, n > 0, are all gradient fields (see, 
for example, [12]). 



Definition 1.5 A linear skew- symmetric operator J : T*(M) — » T(M) is called a Hamil- 
tonian operator or to be Hamiltonian, if for all a, /3, 7 G T*(M), we have 

< a,J'{u)[J(u)(3]-f > +cyde(a,/3,7) = 0. (1.9) 

Its Poisson bracket is defined by 

tiff SH 

{H u H 2 }j{u) =< £_!(«), >, (1.10) 

Su Su 

where H x , H 2 G C°°(M). 4 pair of operators J,M : T*(M) -> T(M) is caZZed a Hamilto- 
nian pair, if cJ + dM is always Hamiltonian for any constants c, d. 



When J : T*(M) -> T{M) is Hamiltonian, we have @ 



J(u) (it) = [J(u)— — (u), J{u)—— (u)\, 



Su 



Su 



Su 
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where Hi,H 2 G C°°(M). This implies that the operator J 4- is a Lie homomorphism 
from the Poisson algebra to the vector field algebra. Moreover if J, M : T*{M) — > T{M) 
is a Hamiltonian pair and J is invertible, then $ = MJ _1 : T(M) — » T(M) defines a 
hereditary operator |l3|, 0]. 

Definition 1.6 ^4 linear skew- symmetric operator : T{M) —* T*(M) is called a sym- 
plectic operator or to be symplectic, if for all K,S,T G T{M), we have 

< K(u),e'{u)[S{u)]T(u) > +cycle{K, S, T) = 0. (1.11) 



If : T(M) — > T*(M) is a symplectic operator, then its second-degree tensor field 
Te G T^M) can be expressed as 

T e = d-y with < 7 (u),if(u) >= f < e{Xu)Xu,K{u) > dX, K G T(M), 

JO 

where d'y is determined by (1.4). It is not difficult to verify that the inverse of a symplectic 
operator is Hamiltonian if it exists and vice versa. We also mention that Hamiltonian and 
symplectic operators can be defined only in terms of Dirac structures E 



Definition 1.7 A system of evolution equations u t = K(u), K G T(M), is called a 
Hamiltonian system or to be Hamiltonian, if there exists a functional H G C°°(M) so that 

<)H 

u t = K{u)=J(u) — (u). (1.12) 
ou 

It is called a bi- Hamiltonian system, if there exist two functionals i?i,i?2 G C°°(M) and 
a Hamiltonian pair J,M : T*(M) — > T(M) so that 

u t = K{u) = J{u)°—±{u) = M(u)°—±{u). (1.13) 
du du 

There is the other kind of Hamiltonian systems, which can be defined by symplectic 
operators. However, the above definition has more advantages in handling symmetries and 
conserved functionals. For a Hamiltonian system Ut = J(u)^(u), the linear operator J-^ 
maps a conserved functional into a symmetry. For a bi-Hamiltonian system, we will be able 
to recursively construct infinitely many commuting symmetries and conserved functionals 



for the system, if either of two Hamiltonian operators is invertible [15] 



In what follows, we would like to develop a theory for constructing integrable cou- 
plings of soliton equations, by analyzing integrable properties of the perturbation systems 
resulted from perturbed soliton equations by multi-scale perturbations. The paper is or- 
ganized as follows. In Section 2, we first establish general explicit structures of hereditary 
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operators, Hamiltonian operators and symplectic operators under the multi-scale pertur- 
bations. We will go on to show that the perturbations preserves complete integrability, 
by establishing various integrable properties of the resulting perturbation systems, such 
as hereditary recursion operator structures, Virasoro symmetry algebras, Lax representa- 
tions, zero curvature representations, Hamiltonian formulations and so on. In Section 3, 
the whole theory will be applied to the KdV equations as illustrative examples. This leads 
to infinitely many integrable couplings of the KdV equations, which include Hamiltonian 
integrable couplings possessing two different hereditary recursion operators and local bi- 
Hamiltonian integrable couplings in 2 + 1 dimensions. Finally, some concluding remarks 
are given in Section 5. 



2 Integrable couplings by perturbations 

2.1 Triangular systems by perturbations 



Let us take a perturbation series for any N > and r > 0: 

N 

un = £%, rji = 7/i(y , 2/1,2/2, ■ • • , y r ,t), Vi = e l x, t e R, x G R p , < i < N, (2.1) 

i=0 

where e is a perturbation parameter and rji, < i < N, are assumed to be column vectors 
of q dimensions as before. When r > 1, ( |2.1|) is really a multi-scale perturbation series. 
We fix a perturbed vector field K = K(e) € T(M) which is required to be analytic with 
respect to e. Let us introduce 



K® = K®(r} N ) ={K(u,e)) {i) (fj N ) = - 



1 &K(u N ,e) 



, < i < N, 



(2.2) 



6 = 



u de % 

where fj^ = (t)q , rj[ , ■ ■ ■ ,rjJf) T as before, and then define the iV-th order perturbation 
vector field on Mn- 

(pev N K)(fj N ) = K n (t)n) = (K^ T (fjN),K^ T (fj N ), ■ ■ ■ ,K^ T (f) N )f. (2.3) 

Here the vector fields on M are viewed as column vectors of q dimensions, and the vector 
fields on Mn, column vectors of q(N + 1) dimensions, as they are normally handled. Since 
we have 

i j 
, ik = ^£ k rik, Uj = ^2e k Vk, < i < j < N, 



d l K( Ui ,e) 



& l K( Uj ,e) 



Bp 1 

it is easy to find that 



e=0 



fc=0 



k=0 



(2.4) 
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Thus the perturbation vector field peic^K = 6 T(M/v) has a specific property that 
the i-th component depends only on 770, 7]i, ■ • • ,r)i, not on any r]j, j > i. 

Let us now consider a system of perturbed evolution equations 



u t = K(u,e), K = K(e) G T(M), 



(2.5) 



where K (e) is assumed to be analytic with respect to e, as an initial system that we start 
from. It is obvious that the following perturbed system 



iiNt = K(un, e) + o{e N ) or UNt = K(un , e) (mod e^), 
leads equivalently to a bigger system of evolution equations 



?7Ar t = K N (fj N ) 



i.e., % = - 



< i < N, 



(2.6) 



(2.7) 



e=0 



where un is defined by ( |2.1| ). Conversely, a solution tjn of the bigger system ( |2.7| ) gives 
rise to an approximate solution un of the initial system (|2.5| ) to a precision o(r^). The 
resulting bigger system (^^) is called an iV-th order perturbation system of the initial 



perturbed system (2.5), and it is a triangular system, owing to (2^). We will analyze its 
integrable properties by exposing structures of other perturbation objects. 



2.2 Symmetry problem 

Let us shed more right on an remarkable relation between the symmetry problem and 
integrable couplings. Assume that a system ut = K(u), K £ T(M), is given. Then its 
linearized system reads as vt = K'(u)[v]. What the symmetry problem requires to do is to 
find vector fields S G T(M) which satisfy this linearized system, i.e., (S(u))t = K'(u)[S(u)] 
when ut = K{u). Therefore (u T , (S(u)) T ) solves the following coupling system 

{ * = * W ' (2.8) 

if S £ T{M) is found to be a symmetry of u± = K(u). This system ( |2.8D has been 
carefully considered upon introducing the perturbation bundle |ffl|. It is the first-order 
standard perturbation system of ut = K(u), introduced in Ref. |5|. Since it keeps complete 
integrability, it provides us with an integrable coupling of the original system m = K(u). 
Therefore the symmetry problem is viewed as a sub-case of general integrable couplings. 

The commutator of the vector fields of the form (K(u), A(u)v) T with A(u) being linear 
has a nice structure: 



K(u) 
A{u)v 



S(u) 
B(u)v 



[K(u),S(u)} 
lA(u),B(u)}v 
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where the commutator [K(u), S(u)] is given by (|1.3|) and the commutator [A(u),-B(u)J of 
two linear operators A(u),B(u) is defined by 

[A(u),B(u)] =A'(u)[S(u)} - B'{u)[K(u)} + A(u)B(u) - B(u)A(u), 

which was used to analyze algebraic structures of Lax operators in fl6|| . Moreover for 
linearized operators, we can have 

lK\u),S'(u)l = T'(u), T = [K,S], K,Se T(M), (2.9) 

which will be shown later on. 



2.3 Candidates for integrable couplings 

Let us illustrate the idea of how to construct candidates for integrable couplings by per- 
turbations. Assume that an unperturbed system is given by 

u t = K(u), K G T(M), (2.10) 

and we want to construct its integrable couplings. To this end, let us choose a simple 
perturbed system 

ut = K{u) + eK(u), (2.11) 

which is analytic with respect to e of course, as an initial system. Obviously this system 
doesn't change integrable properties of the original system ( |2.10D . In practice, we can have 
lots of choices of such perturbed systems. For example, if the system (|2.10| ) has a symmetry 
S G T(M), then we can choose either ut = K(u) + eS(u) or u t = K(u) + e 2 S(u) as another 
initial perturbed system. According to the definition of the perturbation systems in ( |2.7| ), 
the first-order perturbation system of the above perturbed system ( 2.11| ) reads as 



( r] 0t = K(r) ), 

{ (2-12) 
{ 7 llt =K'(r]o)[7 11 }+K{ri ). 

This coupling system is a candidate that we want to construct for getting integrable 
couplings of the original system ut = K{u). In fact, we will verify that the perturbation 
defined by Q2.1| ) preserves complete integrability. Therefore the above coupling system 
( [2.12| ) is an integrable coupling of the original system ut = K(u), provided that ut = 
K(u) itself is integrable. The realization of more integrable couplings, such as local 2 + 
1 dimensional bi-Hamiltonian systems, can be found from an application to the KdV 
hierarchy in the next section. 
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2.4 Structures of perturbation operators 

Rather than working with concrete examples, we would like to establish general structures 
for three kinds of perturbation operators. The following three theorems will show us how 
to construct them explicitly. For the proof of the theorems, we first need to prove a basic 
result about the Gateaux derivative of the perturbation tensor fields. 



Lemma 2.1 Let X = X(e) G TJ(M) be analytic with respect to e and assume that the 
vector field Sn = (Sq , Sf , • • • , Sjf) T € T(Mn), where all sub-vectors S{, < i < N , are 
of the same dimension. Then the following equalities hold: 



d l X(u,N,e) 



E=0 



(w)[s. 



N 



Qi 

de 1 



N 



X'{u N , e) , 0<i<N. (2.13) 

£ = j = Q 



Remark: Note that in ( |2.13| ), we have adopted the notation 

X'(u N ,e)[K(u)] = (X(u N ,e)y(u N )[K(u)}, X = X{e) e TJ(M), K G T(M), (2.14) 

in order to save space. The same notation will be used in the remainder of the paper. 
Proof: Let us first observe Taylor series 

N 



X(u N ,e) = }^- 



i=0 



+ o(e N ). 



6 = 



It follows that 

(X(u N ,s))'(fj N )[S N ] = 
Secondly, we can compute that 
(X(u N ,e))'(rj N )[S N ]= ^ 



N 

E 

i=0 



e i ( WXiuN, e) 



£=0 



(f, N )[S N ]+o(e»). 



6=0 



N N 

X(u N + 5 jS^e) = X'{u N ,e) [j] ePSj 



j=0 



j=0 



A combination of the above two equalities leads to the required equalities in ( 2.13j ), again 
according to Taylor series. The proof is completed. I 



Theorem 2.1 // the operator $ = $(e) : T{M) — > T(M) being analytic with respect to e 
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is hereditary, then the following operator $tv : T{Mj^) — ► T(Mjv) defined by 
(**Kw)) i7 - 



i,j=0,l,- ,N 



1 &- j §(u N ,e) 



(i-j)l dei-i 



5=0. 



$(-%,£) | e=0 



g(AT+l)xg(AT+l) 





J_ d$(u N ,e) 
1! <9e 



e=0 



$({tAr,e)| e=0 



1 ^^(ujy.e) 



AT! de^ 



1 a$(Miv,£) 



e=0 



1! de 



e=Q 



£ = 



(2.15) 



is also hereditary, where un is a perturbation series defined by ( jg. iy 



Proof: Let K N = (K^K?,--- ,i^) T , S n = (Sg\ Sf, • • ■ ,S^) T G T(Mjv), where the 
sub- vectors Ki,S{, < i < iV, are of the same dimension. Since <&n(,Vn) is obviously 
linear, we only need to prove that 



-& n (w)[&n(w)Sn]Kn + ^N(flN)^' N (VN)[SN]K N = 0, 
according to Definition Ol. In what follows, we are going to prove this equality. 



(2.16) 



First, we immediately obtain the i-th element of the vector field ^n(t]n)Kn and the 
element in the (i,j) position of the matrix ^' N (fj^f)[K^]: 



1 tf-i$(u N ,e) 



3=0 



(i-j)l dei-i 



Kj, < i < N, 



e=0 



[& n (J}n)[K n ]).. 
1 d l ~ j 



1 { d l ~i<S>(u N ,e) 



e=0 



(W)[K N ] 



(i-j)l de 1 -! 



ij {i- j)! V d£i ~ j 

N 

&(u N ,e) [J2 e***]), 0<i,j<N, 



e=0 



k=0 



the last equality of which follows from Lemma 2.1 
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Now we can compute the i-th element of & N (fjN)[$N(VN)KN]SN as follows: 
($jv (fyjv) [$N (f) N )K N ] S N )i 

V ' 1 0* '<\Hu x .~ 



E 

i 

E 

3=0 
i 

E 

j=0 

i 

E 





1 


Qi-3 










1 


gi-j 


(* 


-i)! 






i 


Qi-j 


(< 


-i)! 


de*-i 




l 





3=0 

E 



e=0 



e=0 



e=0 



e=0 



k=0 1=0 



(k-l)l de k ~ l 



Ki 



e=0 



1=0 k=l v ' 
N 

&(u N , e) []T e l ($(%, e) + o{e N ~ 1 )) K x 



Ki 



e=0 



s, 



S; 



Sj 



1=0 
N 



1=0 



S; 



1 sh'- 1 



$'{u N ,£) $(u N ,e)Ki 



e=0 



Sj, 0<i<N. 



0<j+l<i 

On the other hand, we can compute the i-th element of ^j^(fj]\f)^' N (fi]\f)[K]\f]Sj^ as follows: 
(® N (f) N )& N (f) N )[K N ]S N ). 



1 &-j$(u N ,e) 



1 &-*$(u N ,e) 



E 

3=0 

i 

j=0 y J ' 

i i j-k 

EEE 

k=0 j=k 1=0 

i i—k i 



E 



1 &>- 



(j - k)\ dei- k 



e=0 



N 



1=0 



■I I l ' | Qj-k-l 



e=0 



1 



-k-l 



(i - - k - l)\ de*-i 

1 &- j $(u N ,E 



e=0 



-k-l 



e=0 



EEE 

k=0 1=0 j=k+l 
i i—k 

^^Ti-k-l)\ de 1 -**- 1 
k=0 i=o v ; 



(i-j)l(j-k-l)\ de*-i 

I Qi—k—l 



-k-l 



e=0 



dei- k - 1 



e=0 



E 



i a 



i—k— l 



0<k+l<i 



(i-k-iy. de*-*- 1 



e=0 



e=0 



&{u N ,e)&{u N ,e)[Ki[Sk 
^(u N ,e)^(u N ,e)[Ki]S k ), 0<i<N. 



Therefore, it follows from the hereditary property of <&(u,e) that each element in the 
left-hand side of ( 2.16| ) is equal to zero, which means that (|2.16| ) is true. The proof is 
completed. I 



Theorem 2.2 If the operator J = J(e) : T*(M) — ► T(M) being analytic with respect to e 



13 



is Hamiltonian, then the following operator Jjy '■ T*{M^) — > T(M^) defined by 
(per N J)(fj N ) = J N (fj N ) 



i.j 0,1,- ,n 



1 



d i+ i~ N J(u N , 



6=0. 



(i+j-N)\ de*+i-*r 

J(un,£)\e=0 
J_ dJ(u? 
1! 9e 



g(JV+l)xg(JV+l) 



£ = 



J(un,e)\ £=0 ^^f^ 



6=0 



J_ d N J(u N ,e) 
W. de N 



£=0 



(2.17) 



is also Hamiltonian, where un is a perturbation series defined by (2.1). 

Proof: Let a N = (a%,aj,--- ,a%) T ,f3 N = ($,(3f,--- ,(3^) T ,j N = (rffinf,--- ,1n) T G 
T*(Mjv), where the sub-vectors o^, 7i, < i < N , are of the same dimension. It suffices 
to prove that 

< a N , J' N (fj N )[J N p N ]j N > +cycle(a N ,(3 N ,^ N ) = 0, (2.18) 

since there is no problem on the linearity and the skew-symmetric property for the per- 
turbation operator Jn(vn)- 

First, based on Lemma |2.l| , we can calculate the element in the (i,j) position of the 
matrix j^(^Ar)[jjv(f7Jv)/^v] as follows: 

N 



1 




Qi+j- 


-N 


(*' + 3 ~ 


N)\ 




-N 


1 




gi+j- 


-N 


(* + 3 ~ 


N)\ 


de*+3 


-V 


1 




Qi+j- 


-N 


(* + 3 ~ 


N)\ 




-V 


1 




Qi+j- 


-N 


(* + 3 ~ 


N)\ 




-V 


1 




Qi+j- 


-N 



i + j -N)\ de i+ i- N 



£ = 



£ = 



£ = 



£ = 



£ = 



j'(u N ,e) (JnPn) 

1=0 

N N 

f(u N ,e)[j2e l 

1=0 k=N- 
N N 

J'(u N ,e)[£e N ~ k ( £ 



d k+l - N J{u N ,e) 



(k + l-N)\ de k+l - N 



Pk 



6=0 



z k+l-N d k+l - N J(u N ,£) 



k=0 

N 



l=N-k 



(k + l-N)\ de k+l ~ N 



6 = 



J'(u N ,e) [J2 e N ~ k {J(uN,e) + o(e k ))(3 k 

k=0 

N 

J'(u N ,e) [j2e N - k J(u N ,e)P, 



k=0 



N 

£ 

k=0 



Qi+j-N 



(i + j-N)\ de i+ i- N 



e 



N—k jl 



J'(u N ,e)[J(u N ,e)p k ]) 



6=0 



N 

k=2N-(i+j) 



Qi+j+k-2N 



+ j + k-2N)\ de i+ i +k - 2N 



j'(u N ,e)[j(u N ,e)p k ], < i,j < N. 



6 = 
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In what follows, let us give the remaining proof for the case of 

Vi = Vi(yo,yi,t) = Vi( x ,£x,t), < i < N. (2.19) 
Suppose that the duality between cotangent vectors and tangent vectors is given by 

<a,K>=[ a T K dx, x £ W, a £ T*(M), K £ T(M). (2.20) 

Jw 

Let us consider the case of x £ II without loss of generality. For brevity, we set 



F 



jk(aN,(3N,7N,d x ) = {aj J\uN,e)[j(u N ,e)f3k\-fj + cycle (ai,/3fc,7j)Y < i,j, k < N, 

(2.21) 



where d x = d yo + ed yi , owing to (^19| ), and we assume that the original Hamiltonian 
operator J(u,e) involves the differential operator d x . Then we can have 



< a N , J' N (f) N )[J N j3 N ]-i N > +cycle(a N ,(3 N ,j N ) 

I Qi+j+k-2N 



oo rco 



oo J — oo 



E 



2N<i+j+k<3N 



(i + j + k - 2N)l d £ i+J+k-2N 



Fijk (atN ,(3n,7n,9 x ) dyodyi 



£ = 



In order to apply the Jacobi identity of J(u,e), we make a dependent variable transfor- 
mation 



Vo = P, yi = Q + £P, 



from which it follows that 



dp — d yo + ed yi , d q — d yi . 
Now we can continue to compute that 

< ajv, JNiVN^JN^N^N > +cycle(a N ,(3 N , ; y N ) 

I Qi+j+k-2N 



(2.22) 



(2.23) 



oo roo 



oo J — oo 



E 



2N<i+j+k<3N 



+ j + k-2N)\ de i+ i+ k - 2N 



e=0 



det 
1 



dyo dyo 
dp dq 

dyi dyi 



dpdq 
Qi+j+k-2N 



(i + j + k-2N)\ de i+ i+ k - 2N 

00 2N<i+j+k<3N v J > 



Fijk(aN,PN,lN,d p ) 

E 

dq = 0. 

In the last but one step, we have utilized the Jacobi identity of J(u,e). 



Fijk(atN ,(3N,lN,d p )dp) dq 



e=0 J -oo 
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The method used here for showing the Jacobi identity can be extended to the other 
cases of the perturbations. Therefore the required result is proved. I 

Similarly, we can show the following structure for the perturbation symplectic opera- 
tors. 



Theorem 2.3 If the operator G = 9(e) : T(M) — » T*(M) being analytic with respect to 
e is symplectic, then the following operator Qn '■ T(M^) —* T*(Mn) defined by 



(per N Q)(f) N ) = Q N (f]N) 
©7v(w)L 



i,j=0,l,™ ,N 



d N - i -iQ(u N ,e) 



(N-i-j)\ de N - i -i 



£=0- 



q(N+l)xq(N+l) 



J_ d N 0(u N ,e) 



1 d N &(u N ,e) 
TT de N 



£=0 



j_ de(ii N ,e) 

1! de 



£ = 



£=0 

@{u N ,e)\ £=0 



Q(uN,e)\ £=0 



(2.24) 



is also symplectic, where un is a perturbation series defined by (2.1). 



2.5 Integrable Properties 

In this sub-section, we study integrable properties of the perturbation systems defined by 
( [2.7D , which include recursion hereditary operators, .ff -symmetries (i.e., time independent 
symmetries), master-symmetries, Lax representations and zero curvature representations, 
Hamiltonian formulations and etc. Simultaneously we establish explicit structures for con- 
structing other perturbation objects such as spectral problems, Hamiltonian functionals, 
and cotangent vector fields. 



Theorem 2.4 Let K = K(e) G T(M) be analytic with respect to e and assume that 
$ = 3>(e) : T(M) — > T(M) is a recursion operator ofu t = K(u,e). Then the operator $jy : 
T(M N ) -» T{M N ) determined by fiudj) is a recursion operator of the perturbation system 
fjjsft = Kn(vn) defined by ( \2. > \ ). Therefore if Ut = K(u,e) has a hereditary recursion 
operator e), then the perturbation system fjm = ^iv(w) has a hereditary recursion 
operator ^>^(fj^). 
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Proof: Let S N = {S% , Sf , ■ ■ ■ , S^) T G T(M N ), where the sub-vectors S i} 0<i<N, are 
of the same dimension. By Lemma Bjl we can compute that 



de k 



e=0 



(W)[K N ] 



Qk 



A' 



Qk 

de k 



e=0 



&(u N ,e)[K(u N ,e)+o(e N )] 



d k $'(u N ,e)[K(u N ,e)} 



de k 



< k < N, 



e=0 



and 



N 



E 

3=0 



£ = 

1 a ; -^'(u W ,£)[S 3 ] 



(i-j)! 9e«-J 



fe=0 

, < % < N. 



e=0 



Therefore, immediately from the first equality above, we obtain the i-th element of <&' N (fjN) [K^]Sn 
as follows: 



{fr N (r) N )[K N ]S N ). = J2 



1 &-> ($'(%, e) [K (u N ,e)]Sj) 



3=0 



, < i < N. (2.25) 



£=0 



Based on the second equality above, we can make the following computation: 



(^(to)[*w(w)Sjv]), 



E 

fc=0 
i 

E 

fe=0 

i i 

EE 

3=0 k=j 

i i 

EE 

3=0 k=j 



1 d l - k 



(i-k)l de { - k 



(K'(u Nl e)[($ N (f) N )S N ) k ]) 



£=0 



i a 



i—k 



(i-k)l 9e*-* 



'(u N ,e) ^2 



1 a fc - J '$(uAT,e) 



(i-Jb)!(fc-j)! de*~ k 



i 



Qi—k 
■i—k 

£ = 

£ fc-3 



S; 



£=0 



(i-j)! s^-j 



£ = 



ik-j) 



K'(u N ,e) 
K'{un, e 



d k - j $(u N ,e) 



de k -i 
d k -i$(u N ,e 



s, 



£ = 



de k ~i 



£=0 



Si 



E 

3=0 

i 

E 

3=0 

i 

E 

3=0 



1 d^' 



(i - j)! de^' 



i a l - j 



K 



£=0 k=j 



{k-j)\ de k ~i 
^ / (n A r,e)[cD(n JV )^+o( e ^)] 



S; 



£ = 



£=0 



1 d i ^{K'(u N ,e)Mu N ,e)S j }) 



(i-i)l 



, < t < iV; 



(2.26) 



£=0 
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[^ N (fj N )K , N (fj N )[S N ]) i 
1 d l - k <S>(u N ,e) 



E 

k=0 



(i - k)\ 



E 



1 d^K'iuN^^Sj 



EE 

3=0 k=j 



1 



9 i - fc $(fi JV ,e) 



ds 



i—k 



de k ~i 
d k ^K'(u N ,e)[Sj 



6 = 



E 

3=0 



1 



i— fc 



< i < N. 



£=0 



e=0 



(2.27) 



e=0 



It follows directly from the above three equalities above fl2.25|) , (|2.26| ) and ( |2.27|) that 
d$ N 



at 



-(Vn)Sn + $ N (r]N)[KN(flN)]SN ~ K N (fjN)[®N(w)SN] + ®N(fjN)K N (fjN)[SN] = 0. 



According to Definition [L^, this implies that the perturbation operator &n{vn) defined 
by fl2.15| ) is a recursion operator of fjNt = ^Ov(^jv)- A combination with Theorem [Oj gives 
rise to the proof of the second required conclusion. The proof is finished. I 

Theorem 2.5 Let K = K(e),S = 5(e) £ T(M) be analytic with respect to e. For two 
perturbation vector fields Kjy, Sn € T(Mjy) defined by (pO), i/iere exists t/ie following 
relation: 

[Kn(vn), S N (fj N )\ = (K N )'(fj N )[S N (fi N )} - (S N )'(fj N )[K N (riN)} = ? N (fj N ), (2.28) 

where T/y € (M^r) is i/ie perturbation vector field of the vector field T(e) = [K(e), S(e)], 
defined by ( \2.3[ ). Furthermore we can have the following: 

(1) if a = a(e) £ T(M) is an n-th order master- symmetry of the perturbed system ut = 
K(u,e), then cjjv £ T(Mjv) defined by \2.3Q is an n-th order master- symmetry of the 
perturbation system fjm = Kn(jin) defined by \2. 



(2) the perturbation system r)m = Kn(vn) defined by (2.7) possesses the same structure 
of symmetry algebras as the original perturbed system u t = K(u,e). 



Proof: As usual, assume that 



T 



S t = (S® T , SW, .., S^F, ft = («, • • • , nf ) 



By the definition of the Gateaux derivative, we first have 



u. 



£>V, 0<i<N. 



k=0 



K{u i + 5^e k S^ k \e) 

5=0 k=0 



d_ 

85 



K(u i + 6S(ui,e) + 8o(£ i ),e) 



5=0 



K'(ui,e)[S(ui,e)] +o(e*), < i < N. 
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Let us apply the equality above to the following Taylor series 

. . £ k d k K(ui,e) , ... 

K (^ s )=J2y d£ k +o(e J ), 0<z<iV, 

k=0 ' £=0 

and then we arrive at 



K' (ui,e)[S(ui,e)} 



k=0 



[Uj,e) 
k\ V de k 



e=0 



fe)[S«] + o(e ,: ), 0<i <N. 



Taking the i-th derivative with respect to e leads to 



(K'( U ,e)[S(u,e)]) (l) (^) =((K(u,e)) w ) 0<i<N. 



(<)Y 



(2.29) 



Now it follows from ( 2.29| ) that for the i-th element of T/v we have 

(T( Uj e))W = (K'(u,e)[S{u,e)])®(rji) - {S'(u,e)[K(u, e))) ( %) 
= ((K(n,e)) (i) )'(^)[5,] - ((5(n,e)) W )'(^)[^] 
= (^(w)!^]), - (S^at)^]),, < i < N. 



This shows that (2.28) holds. All other results is a direct consequence of ( gUp . The proof 
is completed. I 

The relation ( p. 28 ) implies that the perturbation series (2.1) keeps the Lie product of 
vector fields invariant. In particular, the second component of (2.28) yields the Lie product 
property fl2.9|) of linearized operators. In what follows, we will go on to consider Lax 
representations and zero curvature representations for the perturbation system defined by 
( |2.7|) . In our formulation below, we will adopt the following notation for the perturbation 
of a spectral parameter A: 



N 

Atv = £ % fJ-i, AiV = (MO, Ml) - " " > f J 'N) T , 
i=0 



(2.30) 



which is quite similar to the notation for the perturbation of the potential u. Here /ij, < 
i < N, will be taken as the spectral parameters appearing in the perturbation spectral 
problems. A customary symbol \J X ^, x € R, p , will still be used to denote the gradient of 
the spectral parameter A with respect to x. 

Theorem 2.6 Let K = K{e) £ T(M) be analytic with respect to e. Assume that the 
system ut = K(u, e) has an isospectral Lax representation 

j L(u, e)4> = A0, 



4>t = A(u,e)4>, 



(VxA = 0, x e R p ), i.e., {L(u,e)) t = [A{u,e), L(u,e)}, (2.31) 
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where L and A are two s x s matrix differential operators being analytic with respect to u 
and e. Define the perturbation spectral operator Ln and the perturbation Lax operator An 
by 



(per N B)(f) N ) = B N (f) N ) 
B N (f) N )) 



1 1 j 



i,j=0,l,- ,N 



1 d^BiuN^) 



(i-j)l de*-* 



£=0. 



B(u N ,e)\ e=Q 



s(N+l)xs(N+l) 




J_ 9B(njy,e) 
IT We 



e=0 



B(u N ,e)\ e=0 



1 d N B(u N ,e) 



£ = 



1 dB(iiN,e) 
TT Si 



£ = 



B(u N ,e)\ e=Q 



, B = L,A, (2.32) 



where un is given by (2.1). Then under the condition for the spectral operator L that 

if L'(u N )[S(u N )] = o(e N ), S € T(M), then S(u N ) = o(e N ), (2.33) 



the N-th order perturbation system fjwt = ^i\r(?7jv) defined by (2."i) has the following 
isospectral Lax representation 

(L N (rj N )) t = [A N (r} N ),L N (fj N )], (2.34) 

which is the compatibility condition of the following perturbation spectral problem 
( L N (fj N )(f) N = \(f) N , 



.Vj/qA— Vyi^ - 
L wt = A N (r] N )(p N , 

or the following perturbation spectral problem 

( L N (f) N )(f) N = A(/>tv, 



= 0), 



(2.35) 



4>Nt = A N (fj N )(j) N , 



(2.36) 



where the matrix A reads as 

Hah 



A 



L s = diag(l,l,--- , 1) 



(2.37) 



. MiV-fs • • • l^lls Hols 

with the spectral parameter ^j, < i < N, satisfying 

Vy k Hl = 0, < % < N. 



(2.38) 



k+l=i 
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Proof: We first observe that the perturbed system 

u Nt = K(u N ,e) + o(e N ), 



(2.39) 



which engenders precisely the perturbation system fj^t = -Kjv(w) defined by (|2.7| ), Noting 
that L(u, e), A(u, e) are analytic with respect to u and e, it follows from ( |2.31|) that (|2.39| ) 
is equivalent to the following 



Qk 

de k 



e=0 



((L{u N ,e))t - [A(u N ,e),L(u N ,e)]) = 0, < k < N, 



(2.40) 



by use of ( 2.33j ). 

What we want to prove next is that ( 2.40| ) is equivalent to ( p.34| ). Let us compute 
the elements of the differential operator matrix [A^(fj^), ijv(w)]- It is obvious that 
[A^f(fj^f), L^f(fj^f)] is lower triangular, that is to say, 

([A N (fj N ),L N (r) N )})ij = 0, < i < j < N. 
For the other part of [An(t]n), Ljv(f)jv)], we can compute that 

1 d^ k A(u N ,e) 1 d k ~iL(u N ,e) 



(A N (fj N )L N (fj N ))ij = } y 



k=j 



(i-k)\ ds i ~ k 



£ 

k=j 



i — k 



(k-j)\ de k ~3 

d k ~ j L(u N ,e) 



de 



i—k 



e=0 



de k ~i 



e=0 



e=0 



1 d l iA(uN,e)L(uN,£) 



ds l -3 



, < j < i < N, 



e=0 



where the (*) are the binomial coefficients. In the same way, we can obtain 



(L N (fj N )A N (fj N )) 

Therefore we arrive at 

([A N (f) N ),L N (f) N )])ij 



1 d l iL(uN,e)A(uN,e) 

13 = jn% 



de*~i 



< j <i < N. 



£ = 



1 d l -i{A(u N , e),L(u N ,e)] 



, < j < i < N. 



e=0 



Now it is easy to find that ( |2.40| ) is equivalent to ( p. 34 ) . Therefore the perturbation system 
defined by ( |2.7| ) has the Lax representation (2.34). 

Let us now turn to the perturbation spectral problems ( 2.35| ) and ( 2.36| ). Obviously, 
the compatibility condition of the perturbation spectral problem fl2.35| ) is the Lax equation 
( p. 34 ), since the spectral parameter A doesn't vary whatever the spatial variables change. 
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Therefore let us consider the compatibility condition of the perturbation spectral problem 
( |2.36| ), First, we want to prove that 



AA N (fj N ) = A N (ij N )A, 



(2.41) 



if the spectral parameters /Uj, < i < N, satisfy ( p. 38 ). Notice that the condition (|2.3 
on the spectral parameters m, < i < N, is required by 



N 



8=0 1=0 

which is a perturbation version of \J X X = 0. Therefore we have 



A(u N , s, X7 x )\n = X N A(u N , s, yj + o(e N ) 



This guarantees that 
1 &-* 



e=0 



A(un,6, \/ x )\n 



Qi-3 



e=0 



X N A(u N , £, Va,) , < j < i < N, 



which exactly means that the equality ( 2.4l| ) holds. Now we can compute from ( |2.36| ) that 

{L N (f) N )) t lj)N + L N {fi N )A N (fi N )4) N = AA N (fj N )<p N = A N (fi N )A(i) N = A N {fi N )L N {fi N )4> N . 

It follows that the compatibility condition of the perturbation spectral problem ( p.36|) is 
also the Lax equation ( [2.34 ). The proof is completed. I 

The perturbation spectral operator Ln is very similar to the perturbation recursion 
operator $jv, in spite of different orders of matrices. Actually, we may take any recursion 
operator $ as a spectral operator and the system m = K(u) can have a Lax represen- 
tation §t = This Lax representation is usually non-local, because most recursion 
operators are intego-differential. We also remark that two perturbation spectral problems 
above are represented for the same perturbation system defined by ( |2.7| ), which involve 
different conditions on the spectral parameters. For the case of 



N 



N 



u N = ^V77i(x,y,i) = y^£ l rjj(x,£x,t), x € R, 



(2.42) 



i=0 



i=0 



the condition ( 2.38 ) can be reduced to 

Wte = 0, m x + (ii-i >y = 0, 1 < i < N. 



(2.43) 



In the following theorem, a similar result is shown for zero curvature representations of 
the perturbation systems defined by (E 
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(2.44) 



Theorem 2.7 Let K = K(e) G T(M) be analytic with respect to e. Assume that the 
initial system Ut = K(u, e) has an isospectral zero curvature representation 

j 4> x = U(u,X,e)<j), 

{ 4> t = V(u,\,e)(p, 

i.e., (U(u,\,E)) t - (V(u, A, e)) x + [U(u, X,e),V(u, X,e)] = 0, (2.45) 

where U and V are two sxs matrix differential (sometimes multiplication) operators being 
analytic with respect to u, A and e. Define two perturbation matrix differential operators 
Un and Vn by 

(per N W)(fj N ) = Wn(w,P>n) = Wn(t)n) 

1 d'-iWiuN,^)^ 



i,j=0,l,- ,N 



W(f) N ,\ N ,e) 

1 (9Vy(ttjv,Ajy,e) 
TT de 



e=0J s(N+l)xs(N+l) 




e=0 



e=0 



W(f) N , X N ,e) 



e=0 



1 d N W(u N ,X N ,e) 



AT! 



1 dW(u N ,X N ,e) 



e=0 



1! 



de 



e=0 



W(r) N ,\ N ,e) 



e=0 



where W = U,V , and un and \n are given by and ( \2.3(\) . Then under the condition 

for the spectral operator U that 

ifU'{u N )[S(u N )} = o(e N ), S G T{M), then S(u N ) = o{e N ), (2.47) 

the N-th order perturbation system fjm = Kn{t)n) defined by \2. has the following 
isospectral zero curvature representation 



>JV, 



/jr<f>Nvi = u N (fj N ,fi N )(f) I - 

i=0 

4>Nt = v N (fi N ,fi N )(f> N , 

r 

i.e., (U N (VN))t ~ Y, Ui (VN(flN))y t + [UN(fjN),V N (fj N )} = 0, 



i=0 



where the matrix H is defined by 

IsN 



n 



, I sN = diag(I s , ■■■ ,I S ) = diag(l, ■■■ , 1) 
J s(N+l)xs(N+l) N sN 

and the spectral parameters fj,{, < i < N, satisfy 



d Vk m = 0, < i < N. 



(2.48) 



(2.49) 



(2.50) 



(2.51) 



(2.46) 



k+l=i 
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Proof: Note that by use of ( p. 47 ), the zero curvature equation 

((U(u, A, e)) t - (V(u, A, e)) x + [U(u, A, e), V(u, A, e)} = 
for the system ut = K(u, e) yields an equivalent representation 

r 

(U(u N ,X N ,e)) t - y^ y E i (V(u N ,X N ,e)) Vi + U(u N , X N ,e),V(u N , X N , 



i=Q 



= U'(u N )[K(u N ,e)] -^2e l (V(u N ,X N ,e)) yi + [U(u N ,X N ,e),V{u N ,X N ,e)] 



i=0 



o(e N ) (mod e 



(2.52) 



for the perturbation system f\m = -^jv(^at)- in order to recover u^ t = K(un,e) + o(e N ) 
from ( p. 52 ), we need to keep the spectral property X x = d x X = under the perturbation 
up to a precision o(e ). This requires 



N 



8 X X N = o(e N ), d x = ^2e l d yi ,X N = ^eVi> 



i=0 



i=0 



which generates ( |2.51| ). Similar to the proof of Theorem 2.6, differentiating the above 
equation ( p. 52 ) with respect to e up to N times leads to the zero curvature equation 
( |2.49| ), and conversely, we have ( |2.52| ) if ( 2.4S ) holds. Therefore the perturbation system 
f/Nt = Kn^n) has an isospectral zero curvature representation ( 2.49 ). 



The other thing that we need to prove is that the zero curvature equation (2.49) is 
exactly the compatibility condition of the perturbation spectral problem ( |2.48| ). From the 
first system of ( 2.48| ), we have 

r 

y~] ir^jvjfc* = u Nt <})N + u N (p N f 

i=0 

From the second system of ( |2.48| ), we obtain 

4>Nt yi = V Nyt (f> N + V N (f> Nyi , < i < r. 
A combination of the above equalities yields 

r 

y^T?(V Nyi (f>N + Vn0n Vi ) = U Nt 4>N + UnVn<Pn- 

i=0 

On the other hand, we have 

N N 



(2.53) 



y^ / ii l v N $ N y i = y~] v N n l 4> Nyi = VnUn&n, 



(2.54) 



i=0 



i=0 
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by using IIVjv = VatII and the first system of (|2.48| ). It follows from (|2.53| ) and ( 2.541) 



that the zero curvature equation (|2.49| ) is the compatibility condition of the perturbation 



spectral problem ( 2.48| ). The proof is completed. I 



If we consider the specific case of the perturbation defined by ( |2.42| ), then the pertur- 
bation spectral problem and the perturbation zero curvature equation, defined by ( 2.48j ) 
and ( |2.49| ), will be simplified to 

\ ~ " - (2-55) 

I <pNt = V N (r)N,p, N )</) N , 

and 

(U N (rj N )) t - (V N (f, N )) x - U(V N {f] N )) y + [U N {fj N ), V N (fj N )] = 0, (2.56) 

respectively. The involved spectral parameters fa, < i < N, need to satisfy a reduction 
( 2.43| ) of the general condition ( |2.51| ) . 



Theorem 2.8 Let K = K(e) £ T{M) be analytic with respect to e. Assume that the 
initial system ut = K(u,e) possesses a Hamiltonian formulation 

5H 

u t = K(u, e) = J(u, e) — (u, e) , 

ou 

where J : T*(M) — » T(M) is a Hamiltonian operator and H E C°°(M) is a Hamiltonian 
functional. Then the perturbation system fjm = Kn{t)n) defined by ([2.1\) also possesses a 
Hamiltonian formulation 

rim = K N (fj N ) = J N (f) N ) S ( p( : r . (w), (2.57) 

OVN 



where the Hamiltonian operator Jn(vn) is determined by (2.17) and the Hamiltonian 
functional per^H = H^ € C°°(M/v) is defined by 



(per N H){r] N ) = H N (j] N ) = — 



(2.58) 



e=0 



N\ de N 

The corresponding Poisson bracket has the property 

{per N H x ,per N H 2 } j = per N {H u H 2 }j, H U H 2 G C°°(M). (2.59) 



Moreover the perturbation systems fjm = Kn{vn) defined by % ) possesses a multi- 
Hamiltonian formulation 

{> fA x i /. Mper N Hx) f 5{per N H m ) 

rim = K N (rj N ) = Jin{Vn) ft (Vn) = ■■■ = JmNiVN) ~ kVn), 

drjN dr] N 
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if Ut = K(u,e) possesses an analogous multi-Hamiltonian formulation 



5H 

u t = K(u,e) = J\(u,£)-r—{u,e) 

ou 



J m {u,e)^^(u,£). 

ou 



Proof: Assume that 7(e) = f^(e) G T*(M). Let us observe that 



8u 

1 d l (J(u N ,e)j(u N ,e)) 
i\ 



e=0 



3=0 



1 a j -V(-U7V,e) 



j\{i-j)l df-i 



e=0 



, 0<i<N. 



E = 



Thus, noting the structure of Jjy, we can represent the perturbation system as follows 

fjm = K N (fj N ) = J N (fj N )j N (fj N ), (2.60) 
where the cotangent vector field 7/v G T*(Mn) reads as 



IN (v. 



N) 



1 d N 1 T (u N ,e) 



1 dV-t-fiu^e) 



N\ de N 

1 dj T {u N ,e) 



=0 '(iV-l)! de"- 1 



e=0 



1! de 



e=0 



,7 T (uAr,e)| e=0 ) • 



(2.61) 



Let us check whether this cotangent vector field 7/v is a gradient field. If it is gradient, the 
corresponding potential functional has to be the following 



Hn(v 



N 



1 N 1 ^7(An7v,e) 



i=0 

1 d N 



e=0 



N\ de N 



e=o J 



< j(Xu N ,e),u N > dX = — 



N\ de N 



e=0 



The cotangent vector field 77V is indeed a gradient field, because we can show that 

S(per N H) 
JNKVN) = 77 (vn)- 



5fj N 



(2.62) 



According to Definition 1.1 and using Lemma 2.1, for any Si G T(M{rn)) we can compute 
that 



< 



5 1 1 d N H(u N ,e) 



dr/t \Nl de N 
1 8 N 



AT! <9e w 



e=0 



.iV! de^ 

1 



1 d N ~ 



(N-i)\ de N ~ i 



e=0 

H'(u N ,e)[Si(r]i)} 



(N-i)\ de N ~ 



e=0 

1 



5H 



< j(u N ,e),Si(r]i) > 



6 = 



< 



1 d»-+"f(u N ,e) 



(N-i)\ de 



N- 



.Sifa) >, < i < N. 



e=0 
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This equality implies that ( 2.62| ) holds. It follows that the perturbation system ( 2.60 ) is a 
Hamiltonian system. 

Let us now turn to the property ( 2.59 ) for the Poisson bracket. Set 71(e) = ^p-(e), 72(e) 
TTT'(e) € T*(M). In virtue of ( p. 62 ), we can make the computation 



{peT N Hi,peT N H 2 }j N {r] N ) =< {vn),Jn{vn) (Vn) > 



N 

E 

i=0 



N 

E 



< 



d N l ji{u N ,e) 



(N 



in 



8e N ~ 



8fj N 

N 

. E 



Sfj. 



N 



1 



.(i + j-N)\ 



d i+ i- N J(u N ,e) 



g £ i+j-N 



£=0 j=N-i 

1 d N ~^2{uN,e) 



=0 ( n ~jV- 



de N ~i 



> 



E = 



< 



1 



d N ir yi(u N ,e) 



i=0 

1 d N 



(N-i)\ 



e=0 



1 d%J_(u N , s)j 2 (un , e)) 



9e i 



> 



e=o 



< 7i(ujv,e), J{u N ,e)^ 2 {uN,£) >= (per N {Hi, H 2 }j){fi N ) 



e=0 



This shows that the property ( |2.59| ) holds for the Poisson bracket. 

Further, noting the structure of the perturbation Hamiltonian operators, a multi- 
Hamiltonian formulation may readily be established for the perturbation system. Therefore 
the proof is completed. I 



We should realize that two formulas ( 2.58 ) and ( 2.61 ) provide the explicit structures 
for the perturbation Hamiltonian functionals and the perturbation cotangent vector fields. 
The whole theory above can be applied to all soliton hierarchies and thus various inter- 
esting perturbation systems including higher dimensional integrable couplings may be 
presented. In the next section, we will however be only concerned with an application of 
the theory to the KdV soliton hierarchy. 



3 Application to the KdV hierarchy 



Let us consider the case of the KdV hierarchy 

ut = K n = K n (u) = $ = *(u) = d 2 x + 2u x 8- 1 + 4u, n > 0. 



(3.1) 



Except the first linear equation ut = u x , each equation u% = K(u) (n > 1) can be written 
as the following bi-Hamiltonian equation [15| 



u t = K n = J -A = M- 



Su 



Su 



(3.2) 
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The corresponding Hamiltonian pair and Hamiltonian functionals read as 
J = d x , M = M(u) = 81 + 2{d x u + ud x ), 



(3.3) 



H r , 



H n dx, H n 



uf n (Xu) dX, f n = * n u, n > 0, (3.4) 



oo 

m=0' 



H n {u) -- 

where \& = & = d x + Au — 2d~ l u x . Therefore each equation in the KdV hierarchy (|3 
has infinitely many commuting symmetries {K m }^ =n and conserved densities {H m } 

The second equation in the hierarchy ( |3.1[) gives the following KdV equation 

u t = u xxx + 6uu x , (3.5) 
which serves as a well-known model of soliton phenomena. Its many remarkable properties 



were reviewed by Miura [17]. In our discussion, we are concerned only with bi-Hamiltonian 
formulations and consequent symmetries and conserved densities. The bi-Hamiltonian for- 
mulation of the KdV equation ( |3.5| ) can be written down 



u t = J 



6Hi 

6u 



M 



5u 



with two Hamiltonian functionals 



Hq dx 



1 



-u dx, Hi 



Hi dx 



(-uu xx + u 3 )dx. 



(3.6) 



(3.7) 



It has also an isospectral zero curvature representation Ut — V x + [U, V] = with 



U 



-u — A 




V 



-u xx - 2v? + 2\u + 4A 2 



2u-4A 



(3.8) 



where A is a spectral parameter (see [18| for more information). These two properties will 
be used to construct bi-Hamiltonian formulations and zero curvature representations for 
the related perturbation systems. 

In order to apply the general idea of constructing integrable couplings to the KdV 
equations, let us start from the following perturbed equation 



(3.9) 



i=0 



where on are arbitrary constants and the Si are taken from zero function and K n , n > 0, 
so that the series ( |3.9| ) terminates. To obtain integrable couplings of the n-th order KdV 
equation ut = K n , we need to fix So = K n . Various integrable couplings can be generated 



by making the perturbation defined by (2.1). In what follows, we would only like to present 
some illustrative examples. 
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3.1 Standard perturbation systems 



First of all, let us choose the n-th order KdV equation itself as an initial equation: 

u t = K^{u)=K n {u) 

for each n > 1. In this case, the single scale perturbation un = Y^iLo £l Vi{ x ^) leads to a 
type of integrable couplings: 

Vm = K nN (f) N ), N>0, (3.10) 
which are called the standard perturbation systems of ut = K n and have been discussed 



in [19, 20 1 . These systems have the following bi-Hamiltonian formulations [20] 

v- f \ x>n * t S(per N H n ) - 5(yec N H n -.{) 

rjm = K nN (rj N ) = 9 N rj Nx = J N — = M N — , (3.11) 

or)N or] N 

where the Hamiltonian functionals per j^H n , the hereditary recursion operator <3?jv an d the 

Hamiltonian pair {Jn,Mn} are given by 

1 d N H n {u N ) 



per N H n 



N 



N\ de N ' 

$ o(%) 



o 



Jn 



' 



d x 



d. 







(3.12) 



(3.13) 



M N 



M (vo) Mi(?7i 



M (vo) Mi(r?i) 
M N (r) N ) 



,(3.14) 



with 



Mi = Miirn) = 5 i0 dl + 2(d xVi + Vi d x ), = Qifa) = 6 io d? + 2{d x rj i d~ 1 + rji), 0<i<N. 

(3.15) 

Moreover they have infinitely many commuting symmetries {-fC m 7v}m=o anc ^ conserved 
densities {H mN }™ =0 . 

We list the first two standard perturbation systems of the KdV equation (|3.5[) : 

(3.16) 
(3.17) 





— VOxxx 


+ 677o%x, 


mt 


— Vixxx 


+ 6(7/07/1)0;; 


Vot 


— VOxxx 


+ 67/0 7/ 0a: , 




— Vixxx 


+ 6(7/07/1)2;, 


mt 


— f]2xxx 


+ 67/1 7/i x + 6(7/07/2)3; 
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The first-order perturbation system ( 3.16 ) has the following bi-Hamiltonian formulation 



rjit 



a x 

d x 



5fn 



31 + 2770a + 47/0^ 



5(per 1 H ) 



5fji 



(3.18) 



with i)\ = (??0) r /i) T an d the Hamiltonian functionals 



peri#o = / #01 dx, H 01 = 770771, 



1 



perx-ffi = J H u dx, H u = ^(r] rnxx + VoxxVl) + ^VoVi- 
The second-order perturbation system fl3.17|) has the following bi-Hamiltonian formulation 



f)2t = h 



<5(per 2 #i) (5(per 2 # ) 

— Mo 



5fj2 



Sfj 2 



(3.19) 



with the Hamiltonian functionals 



per 2 # = J H Q2 dx, H 02 = 7] r] 2 + -rj 1 , 

per 2 #i = J Hu dx, Hu = ^(vorj2xx + ViVixx + V0xxV2) + 3??o?/i + ^Im- 

Another example that we want to show is the first-order standard perturbation system 
of the fifth order KdV equation ut = if 2 (u): 

f Vot = r]0,5x + WrjoVOxxx + 20r] 0x r] 0xx + 30r]^r] 0x , 

{ mt = Vi,5x + lOrjoxxxVl + WvoVlxxx + ZOr/OxxVix + ZOvoxVixx + SOvoVoxVi + 30r]%r]i x , 

(3.20) 

where ijq^ x and 771,53;, as usual, stand for the fifth order derivatives of 770 and 771 with 
respect to x. It has the following bi-Hamiltonian formulation 



f (5(periF 2 ) ~ 5(peri^i; 



5fn 



5f}i 



(3.21) 



where the Hamiltonian function peri .Hi is given as before and the Hamiltonian function 
periiJ 2 , given as follows 

1 1 20 10 2 5 o 10 40 o 

periii 2 = -r]oxxxxVl + ^VoVixxxx + yW(te??i + —il Q rn xx + -r) 0x r)i + —^Q X 'q 1 x + y 775771. 

3.2 Nonstandard perturbation systems 



Secondly, let us choose a perturbed equation 

u t = K per (u, e) = K n + aeK n , a = const., 
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as an initial equation for each n > 1. This equation can be viewed as 



u t = K^(u,e) 



J 



5(H n + aeH n ) 
5u 



M 



(J + aeJ)— — = (M + aeM) 
ou 



8u 

SHn-l 

5u 



(3.22) 



Therefore the corresponding perturbation systems also have quadruple Hamiltonian for- 
mulations. We focus on the first-order perturbation system under the single scale pertur- 
bation. It has the quadruple Hamiltonian formulation 



fjlt = Ji 
namely, 

m 



(i)<5(periffA" 
Sfji 



£(1) ^(Peri-gw-l) = j(2)<y(perifli 



(2)> 
n , 



M 



(3.23) 



5 " 
d 
9 

3 ad 



5(peiiH n + aH n (r] )) 
Sfji 

5(peiiH n ) 







M 
M Mi 
M 



<5(peri^ n _i + aff„_i(77o)) 



M Q Mi + aM 



^(periff n _i) 



(3.24) 



where the functional periiJ n , periiJ n _i and the operators Mj are defined by ( p. 12 ) and 
( |3.15 ), respectively. Since two Hamiltonian operators and j[ 2 ^ are invertible, we can 
obtain five hereditary recursion operators for the equation ( 3.23j ): 



WV 1 



1 
a 1 



(1)/ f(2)\-l 



1 
-a 1 



$o 
$i - a<E>o $o 



MfVfV 1 



$ 
$i + a<I>o $o 



$ 
$i $ 



where the operators <3?j are defined by ( p. 15 ). These operator structures suggest two classes 
of hereditary recursion operators for the equation ( |3.23| ) 



Po o 
Pi Po 



, *f\p) 



(3.25) 



where P = (Pq,Pi) t with the Pi being arbitrary constants. They are really hereditary 
operators and recursion operators for the equation ( 3.23| ), which can be verified by direct 
computation or by viewing them as the first-order perturbation operators of the initial 
operators Pq + P\£ and Pq<& + PiE&. Therefore the integrable coupling (|3.23|) of the n-th 
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order KdV equation u t = K n (u) possesses two classes of hereditary recursion operators 
denned by ( 3.25 ). These two classes of operators have the property 

A)7o$o 



*f ) (/3)*^( 7 ) = *^C9)*^( 7 ) 



(2)/ 



(3.26) 



A) 7 o$i + (Poll + /3i 7 o)$o /?o 7 o$o 

for any two constant vectors f3 = (Po,f3i) T and 7 = ( 7 o, 7 i) T , which also shows that their 
product can not constitute completely new recursion operators. 

We can also start from the perturbed KdV type equation 

ut = K per (u,e) = K n + ae j K i3 , a = const., a ^ 0, (3.27) 

where ij is a natural number. Let us illustrate the idea of construction by the following 
specific example 



u t = KP ei (u,e) =K n + ae 2 K n+l (n > 1), 
which can be viewed tri-Hamiltonian system: 



ut = K pcr = J 



5(H n + ae 2 H n+ i) 
Su 



(J + ae 2 M) ^ 

ou 



M 



6(g w -i + ae 2 H n 
8u 



(3.28) 



(3.29) 



Therefore, according to Theorem 2.8, the second-order perturbation system of the the 
perturbed system ( |3.28|) 



r] 0t = K n (j] Q ), 
mt = K' n (rio)[r]i}, 

1 d 2 K n (v,2) 



(3.30) 



e=0 



+ aK n+ i (jiq 



possesses the following tri-Hamiltonian formulation 



V2t = «/; 



(1) SHn 



(1) 



J, 



(2) SH^ 



(2) 



Sfj 2 5f] 2 
with a triple of Hamiltonian operators 



J. 



(3.31) 



J, 



(i) 



d x d x 

d x , J 2 (2) = d x 
<9x J |_ d x qM 

and the corresponding three Hamiltonian functionals 

r(l) 



, J. 



(3) 



M 
M Mi 
M Mi M 2 



(3.32) 



Hn (m) = (per 2 iT n )(f) 2 ) + a-f/n+ifao), 
M 2) fe) = (per 2j ff n )(r? 2 ), 
Hk >(f)2) = (per 2 i? n _i)(r/ 2 ) + aH n (r] ). 



(3.33) 
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Similarly, the perturbation system ( 3.30| ) has also two classes of hereditary recursion op- 
erators: 



A> o o 

01 0o 

02 Pi 00 



® [ 2 2 \0) 



0o^o 

/3 $i + /3i$ /3 $ 

/3 $2 + /9l^l+/32^0 /3 $i + /?i$o 0o$o 

(3.34) 



where the operators $j are defined by ( 3.15| ) and f3 = (0o, 0i, 02) T is a constant vector. 

Let us fix n = 1 and then the system ( |3.30| ) gives an integrable coupling of the KdV 
equation (|3.5| ), which possesses the following tri-Hamiltonian formulation 



V2t 



J, 



(2) 



J, 



5r]2 8m 5fj2 

with three Hamiltonian functionals 

Hi\fi 2 ) = ^(rjoV2xx + ViVixx + Voxxm) + SiloVi + ^VoV2 
+a(\r] Q r]Q XXXX + ™rj%r]oxx + |r?o^ox + T^o)' 
H[ 2 \m) = mm + Ivf + a{^VoVoxx + r/o), 

Hxim) = ^{fioV2xx + mvixx + voxxm) + Zmvi + ^vlm- 



(3.35) 



(3.36) 



In order to distinguish the standard perturbation systems defined by (3.10), the inte- 
grable couplings of the n-th order KdV equation u t = K n , defined by ( 3.23 ) and ( p. 30 ), are 
called the non-standard perturbation systems. Interestingly, each of these systems has both 
a local multi-Hamiltonian formulation and two classes of hereditary recursion operators. 



3.3 2+1 dimensional integrable couplings 



Thirdly, let us consider a case of bi-scale perturbations ( |2.42 ), i.e. 



N 



UN 



^V?7i, rn = r]i(x,y,t), y = ex. 



i=0 



In order to present explicit results for integrable couplings, we take the KdV equation 
(|3.5D as an illustrative example, due to its simplicity. We recall that the KdV equation 
( [3.5D has the bi- Hamiltonian formulation (|3.6| ) and the Lax pair ( |3.8| ), 

Let us introduce the bi-scale perturbation series above into the KdV equation ( |3.5| ) 
and equate powers of e. As a iV-th order approximation, we obtain a 2 + 1 dimensional 
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perturbation systems of evolution equations 



Vot! = Voxxx + SiJoVox, 
Vih = Vixxx + 3rj 0xxy + 6(rj r]i) x + 6r/ r?oj/, 
V2ti = V2xxx + 37?!^ + 3r) 0xyy + 6(r)oV2)x + SviVix + G(voVi)y, 

Vjti — f]jxxx ~\~ ^Vj — l,xxy ~i~ 3T]j—2,xyy ~i~ Vj—3,yyy 

+6( ELo ViVj-i,x + Ei=0 mVj-i-l,y) , 3 < j < N. 



(3.37) 



This system has been already presented in [^(J. It follows from our general theory that it 
gives an integrable coupling of the KdV equation (|3.5[) . 



In what follows, we would like to propose a bi-Hamiltonian formulation and the con- 
sequent hereditary recursion operator for the system ( 3.37| ). To the end, we first need to 
compute a perturbation Hamiltonian pair by Theorem |2.2|: 



Jn 



dx 

d x d y 
d x d y 



M 



N 



dx d y 




(3.38) 



P(s) 



e=0 



P(s) 

1_ dP{e) 
1! de 



le=0 



£ = 



Pie) 



e=0 



1_ dP(e) 
1! de 



e=0 



J_ d 2 P(e) 
2! de 2 



e=0 



P(e) 

1 dP{e) 
1! de 

J_ d 2 P(e) 
2! de 2 



J_ 8 N P(e) 
N\ de N 



6=0 



£=0 



e=0 



(3.39) 



where the differential operator P{e) represents 



P{e) = (8 X + edyf + 2[(d x + ed y )u N + u N {3 x + ed y )}. 

The explicit expressions for various derivatives of P(e) with respect to e can be obtained 
as follows: 

' P(e)| e=0 = dl + 2(d xVo + rtod x ), 

3d%dy + 2{d x r]i + rjid x ) + 2(d y r] + f] d y ), 

= 3d x d% + 2(d x 7] 2 + r) 2 dx) + 2(d y r] 1 + r^dy), (3 40) 

= djj + 2(5^3 + r] 3 d x ) + 2(<9 y ?72 + m^y), 
-■ 2(d x rji + Vidx) + 2(d y r)i-i + m-idy), 4 < i < N, 



1 


dP(e) 




1! 


de 


e=0 


1 


d 2 P(e) 




2! 


de 2 


6=0 


1 


d 3 P(e) 




3! 


de s 


e=0 


1 






t! 


de i 


6 = 
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which gives rise to an explicit expression for the Hamiltonian operator Mn- Secondly, 
we need to compute the Hamiltonian functionals for the system ( |3 .37 ). Note that d x — > 
d x + ed y , and thus, under the perturbation ( p. 42 ), we have 



N 



i=0 



Further, by Theorem 2.5, we obtain two perturbation Hamiltonian functionals: 



perjv#o 
per N Hi 



1 8 N H (u N ) 



N\ de N 
1 d N H 1 (u N ) 



e=0 



dxdy 



1 N 



rjir] N ^i dxdy, 



(3.41) 



£ = 



dxdy 



i=0 



[5 E 



i+j=N 



N\ de N 

+ ^2 ViVjxy + \ ^2 ViVjyy + WiVk] dxd V- ( 3 - 42 ) 



i+j=N-2 



i+j+k=N 



Now the following bi-Hamiltonian formulation for the system Q3.37) becomes clear: 



f}Nt = Jn 



5(per N H 1 



6fj 



M N 



S(per N H ) 



Si) 



N 



(3.43) 



where Jn, Mn, per/vPo and per n H\ are defined by ( 3.38| ), ( |3.39| ), ( |3.41| ) and ( |3.42| ), re- 
spectively. It should be realized that the 2 + 1 dimensional bi-Hamiltonian system ( |3.43| ) 
is local, because the Hamiltonian pair {Jn,Mn} involves only the differential operators 
d x and d v . 



Theorem 2.4 guarantees the existence of a hereditary recursion operator for the system 



( |3.37| ). It is of interest to get its explicit expression. Note that the first Hamiltonian 
operator Jn has an invertible operator 



(Jn)- 1 



Pv Pn-i 
Pn-i 



(3.44) 



Po 
where the operators Pi are defined by 

p = a-\ a = -d~ 2 d y , • • • , p = i-iyd-^di, ■ ■ ■ , pv = (-i)^-"- 1 ^. (3.45) 

Therefore, the corresponding hereditary recursion operator is determined by &n = MnJn , 
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but it can also be computed directly by Theorem 2.1 : 

<f>(u N )\e=0 

<5>(un)\s=o 



N 



1 d<S>(u N ) 
1! de 



e=0 



1 d N <S>(u N ) 
W. W 



e=0 



1 d<S>{u N ) 
TT de 



e=0 



$(%)|e=0 



(3.46) 



Here the operator 3>(&jv) is defined by 

$C%) = (d x + tdy) 2 + 2(uj Va; + eu Ny )(d x + e^) -1 + M N , 
and thus its iV + 1 derivatives with respect to e are found to be 
' $(%)| £=0 = Ql + 2r ?0 x9" 1 + 4r? , 

2<9 a; 3 2/ + 2(7?i x + vo y )d~ 1 - 2rj 0x d x 2 d y + 4t?i, 



1 8<S>(u N ) 
1! 9e 

1 d 2 <3>(u N ) 

1 d k $(u N ) 
k\ de k 



e=0 



£=0 



d 2 + 2(7/2* + 1 - %(vix + TRH/)Sr 2 9 y + 2r] 0x d x 3 d 2 + Arj 2 , 



£=0 



£ i+J=fc (-iy (2r/ M + ,,, .,,)<),■ ' X iK + 4%, 3 < k < N, 



(3.47) 



where we accept r/_i = 0. 

Let us now show the corresponding zero curvature representation for the 2 + 1 dimen- 



sional perturbation system (3.37). By Theorem 2.7 or (2.56), the zero curvature represen- 
tation for the system fl3.37| ) can be given by 

Urn - V Nx - UV Ny + [U N , V N ] = 0, (3.48) 

where three matrices II, Un and Vn read as 




n 



/27V 



, I 2N = diag(/ 2 , ■•• ,h) = diag(l, • • • , 1), (3.49) 

2(7V+l)x2(/V+l) " N ' "^N^ 



U N 







L u N ■■■ U\ u 

with the Ui, V% being determined by 



Vn 



V 

V N ••• 







V x V 



(3.50) 



i\ cV 



V=- 



l#V(u N ,\ N ) 



£=0 



£=0 



-rji - m 
S i0 

2rji - -r] ix - m-l,y 

Qi = -Vixx ~ 2Vi-l,xy ~ Vi-2,yy ~ 2 (r] k T]i - fJb k T]i - 2/J fc /i/), < i < N, (3.53) 

k+l=i 



, 0<i<N, (3.51) 
< i < N, (3.52) 
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where we accept that = t?_ 2 = 0, and U, V are defined by (|3.8|) . Of course, we 
require the condition (|2.43| ) on the involved spectral parameters //j, < i < iV, in order 
to guarantee the equivalence between the system ( |3.37 ) and the zero curvature equation 

(H)- 

In particular, the first-order bi-scale perturbation system 

J Vot = Voxxx + QrjoVox, 
{ Vit = mxxx + Zrioxxy + 6(r/ f?i)x + QrjoVOy, 
has a local 2 + 1 dimensional bi-Hamiltonian formulation 

5 (peri Hi) 



(3.54) 



Sfjx 



y 5(periF ) „ 

Mi — , m 

or/i 



no 
m 








d x 




Jl = 






, Mi = 


_ d x 


Oy _ 



dl + 2 V0x + 4 m d x 

dl + 2t] Qx + 4r] d x 3dld y + 2rji x + 2r] 0y + 4r]id x + 4%^ 



peri#o 



rjorjidxdy, perii/i 



[^rjoVixx + iJoVoxy + t^ViVoxx + ^Ivi) dxdy. 



. (3.58) 



Here the extended variables rjo(x,y,t) and rji(x,y,t) are taken as a potential vector fji. 
Moreover the above Hamiltonian pair yields a hereditary recursion operator in 2 + 1 
dimensions 

d 2 x + 2ri Qx d~ l + 4% 
2d x d y - 2 V0x d- 2 d y + 2(rj lx + r/o^d" 1 + 4r?i d 2 x + 2T/0X9- 1 + 4% 

The system ( p. 54 ) was furnished in p0[ , its Painleve property and zero curvature rep- 
resentation were discussed by Sakovich [21|, and its localized soliton-like solutions were 
found in [f22[|. All these properties show that the system ( 3.54| ) is a good example of typical 
soliton equations in 2 + 1 dimensions. 



(3.55) 

(3.56) 
(3.57) 



4 Concluding remarks 

We have developed a theory for constructing integrable couplings of soliton equations by 
perturbations. The symmetry problem is viewed as a special case of integrable couplings. 
The general structures of hereditary recursion operators, Hamiltonian operators, symplec- 
tic operators, Hamiltonian formulations etc. have been established under the multi-scale 
perturbations. The perturbation systems have richer structures of Lax representations and 
zero curvature representations than the original systems. For example, in the higher di- 
mensional cases, the involved spectral parameters /Xj, < i < N, may vary with respect 
to the spatial variables, but they need to satisfy some conditions, for example, 

= 0, fJL ix + tH-i, y = 0, 1 < i < N, 
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in the 2 + 1 dimensional case of the perturbation 

N N 



^e l rji{x,y,t) = ^e l rn(x,ex,t), x G R. 



u n = 

i=0 i=0 

The resulting theory has been applied to the KdV soliton hierarchy and thus various 
integrable couplings are presented for each soliton equation in the KdV hierarchy. The ob- 
tained integrable couplings of the original KdV equations have infinitely many commuting 
symmetries and conserved densities. Linear combinations of the KdV hierarchy containing 
a small perturbation parameter may yield much more interesting integrable couplings. For 
example, the KdV type systems of soliton equations possessing both multi-Hamiltonian 
formulations and two classes of hereditary recursion operators have been presented and 
what's more, local 2 + 1 dimensional bi-Hamiltonian systems of the KdV type with hered- 
itary structures have also been constructed. 

Our success in extending the standard perturbation cases to the non-standard cases 
and the higher dimensional cases are based on the following two simple ideas. First, we 
chose the perturbed systems as initial systems to generate integrable couplings for given 
integrable systems. The method of construction is similar to that in ||. Only a slight dif- 
ference is that new initial systems themselves involve a small perturbation parameter, but 
importantly, such initial perturbed systems take effect in getting new integrable couplings. 
In particular, our result showed that the following non-standard perturbation system 

Ut = K(u), , , 

v t = K'(u)[v]+K(u), l4 ' ij 

keeps complete integrability. Therefore, this also provides us with an extension of integrable 
systems. Secondly, we took the multi-scale perturbations, by which higher dimensional 
integrable couplings can be presented. Indeed, the multi-scale perturbations enlarge the 
spatial dimensions and keeps complete integrability of the system under study. A concrete 
example of integrable couplings resulted from the multi-scale perturbations is the following 
system 

J VOt = VOxxx + ^VoVox, 
I mt = mxxx + 3ri xxy + 6(770771)2 + 6r] V0y, 
which has been proved to be a local bi-Hamiltonian system. 



A kind of reduction of the standard multi-scale perturbations defined by (2.1) may be 
taken, which can be generally represented as 

N 

un = /j^Vj, Vi = ^(2/0,2/1,2/2, ■ • • ,Vr,t), yj = e l ix, t £ R, x £ R p , < i < N, 
3=0 
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where the can be any two finite sets of natural numbers. This kind of perturbations 



can be generated from the standard perturbations (2.1), if some dependent variables rji 
are chosen to be zero and the other dependent variables are assumed to be independent 
of some dependent variables j/j. They yield more specific integrable couplings. There is 
another interesting problem related integrable couplings. Could one reduce the spatial 
dimensions of a given integrable system while formulating integrable couplings? If the 
answer is yes, it is of interest to find some ways to construct such kind of integrable 
couplings, i.e., to hunt for the second part S(u,v) with v being less dimensional than u to 
constitute integrable systems with the original system Ut = K{u). 

There exist some important works to deal with asymptotic analysis and asymptotic 



integrability p3| , £4|, 25, to which the study of the perturbation systems may be helpful. 
It is also worthy mentioning that our 2 + 1 dimensional hereditary recursion operators, 
for example, the operators denned by ( 3.46| ) and ( p. 47 ), are of the form described only 



by independent variables involved. Thus they are a supplement to a class of recursion 



operators in 2 + 1 dimensions discussed by Zakharov and Konopelchenko 27], and a class 
of the extended recursion operators in 2 + 1 dimensions including additional independent 
variables, introduced by Santini and Fokas ]2!| and Fokas and Santini ||29|| . The other 
properties such as Backlund transformations, bilinear forms and soliton solutions might 
be found for the resulting perturbation systems. A remarkable Miura transformation [ 30 1 



might also be introduced for the 2 + 1 dimensional perturbation systems ( 3.37) ), which will 



lead to new 2 + 1 dimensional integrable systems of the MKdV type. All these problems 
will be analyzed in a further publication. 
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